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Abstract 

Let Ho = [—1,1], and define the segments E n recursively in the 
following manner: for every n = 0, 1, ... , let H n+ i = E n D [a n +l ~ 
1, a n+ i + 1], where the point a n+ \ is chosen randomly on the segment 
S n with uniform distribution. For the radius p n of S n we prove that 
n(p n — 1/2) converges in distribution to an exponential law, and we 
show that the centre of the limiting unit interval has arcsine distribu- 
tion. 

Keywords: Arcsine law; Continuous state space Markov chain; 
Poisson-Dirichlet law; Intersection of convex discs. 

1 Introduction 

We consider the following stochastic process. Let S = [—1,1], and define 
the segments H n recursively in the following manner: for every n — 0,1, ... , 
let 

^n+l = H [a n+ i — 1, a n+ i + 1], 

where the point a n+ i is chosen randomly on the segment S n with uniform 
distribution. After n steps, one obtains the segment 

— [Z n ~ Pn, Z n + p n \. 

The (centre, radius) process (Z n , p n ) is a continuous state space Markov 
chain. The radius sequence (p n ) is monotonically decreasing, and it is easy 
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to see that with probability 1, p n — > 1/2. Moreover, (Z n ) is a convergent 
sequence, assuming values on [—1/2, 1/2]. Denote lim^oo Z n by Z. 
We are interested in the most straightforward questions: 

(1) What is the distribution of Z n and p n for a given n ? 

(2) What is the asymptotic behaviour of the radius? 

(3) What is the limit distribution of the centre? 

Our work was motivated by the following problem formulated by Balint 
Toth (Toth, 2010) some 20 years ago with K being the unit disc of the 
plane. Let K = Kq be a convex body in the M. d that contains the origin, and 
define the process (K n ,p n ), n ^ 1, similarly to the above construction: let 
Pn+i be a uniform random point in K n , and set K n+X = K n fl {p n+ \ + K). 
Clearly, (K n ) is a nested sequence of convex bodies which converge to a non- 
empty limit object, again a convex body in M. d . What can we say about the 
distribution of this limit body? When K is the unit disc the limit object 
is almost surely a convex disc of constant width 1. The present note deals 
with the 1-dimensional analogue of this problem. Intriguingly, apart from 
almost trivial results, nothing is known about the related questions even in 
the plane. 

Another direction of generalising the problem treated here is to choose 
the subsequent centres according to a previously fixed distribution instead of 
the uniform one. Research on this version is currently in progress. 

The paper is organised as follows. In Section 2 we give a recursion for the 
density function of p n , which allows us to explicitly calculate the expectation 
for small n's. In Section 3 we show that n(p n — 1/2) converges in distribution 
to an exponential law, which actually shows the rapidness of the process. 
Section 4 contains the limit distribution of the centre, while in the last section 
we derive a somewhat unexpected connection between the process and the 
Poisson-Dirichlet distribution. 

2 First observations 

At the (n + l)st step, we separate two cases according to the location of 
a n+1 . First, if a n+ i G [Z n + p n — 1, Z n — p n + 1], then no change occurs to the 
segment, and thus S n+ i = S n . In the second case, when a n +i is close to one 
of the endpoints of S n , the centre moves and the length decreases. Introduce 
yet two new random processes measuring the change of the location of the 
centre by 

£nA„ = Z n — Z n _i, n ^ 1, 



2 



with e n = ±1 and X n ^ (if Z n = Z n _i, then of no consequence, let e n = 1). 
Thus, for n ^ 1, 



moreover, 



JO W.prob. (1 - Pn-l)/Pn-l 

} uniform on [0,p„_i - 1/2] w.prob. (2/?„_i - l)/p n _i, 



(e n = l|X n ^0) = i. 



By definition, 

n 

Z n = £i-Xt) (2) 

i=l 

and it is easy to see that 

p n = 1 - ^Jjfj =: 1 - £> n . 



1=1 

\oo 



Thus, with probability 1, ^J°X; = 1/2. By an inductive argument, it 
follows that S n has a continuous distribution. Denote by f n (x) the probability 
density function of S n . Using the Markov property and flTJ it is easy to 
express f n+ \ in terms of f n : for ^ x ^ 1/2, 

nS n+ i >x)= [^nS^ > x\S n = y)f n (y) dy 



whence differentiating 



1-y 



Ui(x) = fn(x) T ^- + 2 r^-dy. (3) 



1 - x Jo 1 - y 
The first few examples are, for ^ x ^ 1/2, 

fi(x)=2 

2x 

f2(x) = 41n(l -x) 

1 — x 

, , , 2x(2-x) 4(1 -2x) , M . An M ,, 2 
(1 — x) z 1 — X 

In order to calculate the expectation of 5 n (and thus Ep n ), we consider 
the Taylor series expansion of /„ about 0: 

oo 

fn(x) ^ C n ^X . 
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Based on ([3]), one readily obtains the formula 



C n +l,k 



k + 2 



fc-i 



i=0 



(4) 



and hence 



ES„ 



1/2 



2 -(fe+2) J 



T^C n _ij ^ -jT> ( 5 ) 



k=0 j=0 k=j+l 



which may be useful for obtaining a recursive expression involving only the 
expectations. These formulas also enable us to efficiently compute the ex- 
pectations for relatively small n's, see the figure below. 
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Figure 1: The first 70 values of n(l/2 - ES n 



3 Asymptotics of the length 

We determine the asymptotic behavior of p n , the radius of the nth segment. 
Theorem 1. We have the distributional convergence 

n { Pn ~l) ^ Ex P( 4 )- 
Accordingly, for any k ^ 1, 
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Proof. Since p n = 1 — S n , it is equivalent to prove the corresponding limit 
theorems for 1/2 — S n . 
It follows from (JT]) that 



s, 



n+1 



W.prob. 5 n /(l - Sn) 

uniform on [£„, 1/2] w.prob. (1 — 2S n )/(l — 



Observe that given S n , S n+ i has the same distribution as ma.x{S n , U}, where 
U is a uniform random variable on [S n /2, 1/2], independent of 5* n . Let 
a ^ 1/2, and C/f, U", ... be independent, uniform random variables on [1/4— 
a/2,1/2]. Define M° = maxi^ n ?7". It is well-known in extreme value 
theory (Billingsley, 1995, p. 192), and easy to check that 

nQ-M") A Exp(4/(1 + 2a)). (6) 

We say that X stochastically dominates (or just dominates) Y if P(X > 
x) ^ F(Y > x) for all x eR. 

To obtain a lower estimate, we use that S n ^ 1/2, and thus an easy 
induction argument shows that M° dominates S n . From this and (jHJ) we 
obtain 

P (n Q - SU > arj > P U Q - M° ] > x 1 ->• e" 4a; 

i.e. 

liminfP ( n f - - S n ) > x ] ^ e" 4x , x ^ 0. 

n->oo \ \2 J J 

The almost sure convergence S n — > 1/2 heuristically means that with 
overwhelming probability, 1/2 — S n /2 « 1/4, hence for n sufficiently large, 
5* n behaves approximately like M° does. This heuristic idea is made precise 
as follows. Fix the small positive numbers > and e > 0. Since S n — >■ 1/2 
a.s., for n sufficiently large, P(S^ n < 1/2 — e) ^ e. Moreover, if ^ 1/2 — e, 
then 5* n is minored by M^_^ n . This can be shown again by induction. Thus, 
by®, 

P ( // ( I - S n ) > x ) = P [ S„ < ~ - - 



2 'V I \ 2 n 

P P ( P., < J - - 

2 n 



< p I < \ ~ I 



4x(l-/3) 

-> e i+ 2 <= + £ 
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i.e. for x ^ 

(/ 1 \ \ 4:r(l-/3) 

n I - — S n J > x I ^ e !+ 2e + £. 

Since e and /3 are arbitrary, this gives the distributional convergence. 

To prove the convergence of the moments, it is enough to show that for any 
k ^ 1, the sequence {n k f | — SVi) is uniformly integrable (Billingsley, 

1 /2 

1995, Theorem 25.12). The fact that S n dominates M n readily implies 



1 I - S„ ) >x) ^e" 2 ' 



2 

which shows uniform integrability. □ 

In particular, Ep n = 1/2 — l/(4n) + o(n~ 1 ) and Var(p n ) ~ l/(16n 2 ). The 
question of obtaining an exact or asymptotic formula for KS n using only the 
recursive relations (j4j) and ([5]) remains open. 



4 Limit distribution 

In this section, we determine the limit behaviour of Z n . Let F(x) denote the 
cumulative distribution function of Z; based on the definition of Z, it follows 
that F(x) = if x ^ -1/2 and F(x) = 1 if x ^ 1/2. 

Theorem 2. The distribution of Z is a translated arcsine law: for —1/2 ^ 
x ^ 1/2, 

2 

F(x) = — arcsin \j x + 1/2 . 

7T 

Proof. By (J2]), we have to determine the limit distribution of ^2 £jXj. Clearly, 
this is not affected by the steps where Xj = 0. We introduce the thinned 
process (Z n , p n ) as follows: for n ^ 1, let £ n be independent Bernoulli (1/2) 
random variables, and X n be a uniform random variable on [0, 1/2— ^Q] ■ 
The centre of the segment after the nth step of the thinned process is given 
by Z n = and the radius is p n = 1 — Y^x^i- Plainly, 

oo 

1=1 

Introduce r n = 2p n — 1. If Z7i, ?7 2 , . . . are i.i.d. Uniform(0, 1) random vari- 
ables, then setting (Z ,r ) = (0, 1), 
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Z n +i = Z n + — £ n +i(l — U n+ i)r n , 
r n+i — U n +ir n . 

Notice that after choosing (Z\, rx), the process (Z 2 , r 2 ), (Z 3 , r 3 ), . . . is a scaled 
and translated copy of the original one, which implies the distributional equa- 
tion 

z = n z 1 + z 1? (8) 

where Z' is independent from (Zi,ri), and has the same distribution as Z. 
Thus, for every x G [—1/2, 1/2], 



l-2yy y VI -2y 
/ ,/ ( n\, <*z + (1 + 2s) H -JM- dz 





- 9xW - >L , , , , 

^ .1-2^)2 v ,y B (i + 2 ^2 

= (1-2*) f F{Z) dz+(l + 2x) T F ^ dz+i±^ 

This also shows that F is continuously different iable, and by differentiating 
we arrive at 

n r F ( z ) j o /" 1/2 ^C*) j 4xF(a;) 1 
F'(x) = -2 / dz + 2 dz H + - . 

Once again, we derive that F is twice different iable (being the reason for 
starting with the distribution function rather than the density function), 
whence 

Taking into account that F'(x) is a density function on [—1/2,1/2] yields 
that the solution is 

F'(x) ' 



7iV(l/2 + x)(l/2-x) ' 
the desired density function. □ 



5 Further remarks 

Setting Vi = UiU 2 ...Ui(l - U i+1 ), implies that Z n = 1/2 J2q~ v^+i, 
hence the limit Z has the infinite series representation 

^ 00 

Z=-J2<i+i- (9) 

i=0 
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It is easy to check that v j = 1 almost surely, thus v := (v , Vi, . . .) G A, 
where 

A = |x = (x ,xi, . . .) : ^2 x i = > 0,i = 0,1,...J 

is the infinite dimensional simplex. The construction of the random vector 
v implies that it has the so-called GEM (Griffiths-Engen-McCloskey) distri- 
bution with parameter 1 (see the residual allocation model in Bertoin, 2006, 
p. 89). This distribution appears in various contexts, such as prime factori- 
sation of a random integer (Hirth, 1997); and in particular, the decreasing 
reordering of the GEM distribution is the so-called Poisson-Dirichlet distri- 
bution, which is one of the most important distributions in fragmentation 
theory, see Bertoin, 2006. 

Using this terminology and OH]), Theorem [5] can be reformulated as: if 
v = (vq,vx, . . .) is a GEM distributed random vector, and £,£1,^2, ... is an 
iid sequence of Bernoulli random variables such that P(£ = 1) = 1/2 = 
P(£ = —1), which is independent of v, then Yl^o Vi ^ i + 1 nas arcsrne distribu- 
tion. This theorem was first proved by Donnelly and Tavare (1987), using 
the construction of the Poisson-Dirichlet distribution by means of an inho- 
mogeneous Poisson process. Later Hirth (1997) also gave a proof by using 
the method of moments. As far as we know, our proof, solving an integral 
equation based on the distributional equality jS]), is new. 
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